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We present a generic formulation of the continuum elasticity of an isotropic crosslinked active gel. 
The gel is described by a two-component model consisting of an elastic network coupled frictionally 
to a permeating fluid. Activity is induced by active crosslinkers that undergo an ATP-activated 
cycle and transmit forces to the network. The on/off dynamics of the active crosslinkers is described 
via rate equations for unbound and bound motors. For large activity motors yield a contractile 
instability of the network. At smaller values of activity, the on/off motor dynamics provides an 
effective inertial drag on the network that opposes elastic restoring forces, resulting in spontaneous 
oscillations. Our work provides a continuum formulation that unifies earlier microscopic models 
of oscillations in muscle sarcomers and a generic framework for the description of the large scale 
properties of isotropic active solids. 



I. INTRODUCTION 

Much recent theoretical effort has focused on mod- 
eling the effect of motor activity on the cell cytoskele- 
ton. The cytoskeleton is a highly heterogenous polymer 
gel, mainly composed of filamentary actin crosslinked 
by a myriad of globular proteins pQ. These include 
proteins that preserve the isotropic nature of the net- 
work (e.g., filamin), proteins that induce bundle for- 
mation (e.g., fascin or vilin), and molecular motor pro- 
teins, such as kynesins and myosins, that are capable of 
transforming chemical energy into mechanical work [2J. 
Motor proteins hydrolyze adenosine-tri-phosphate (ATP) 
and convert it to adenosine-di-phosphate(ADP) and in- 
organic phosphate(P). The free energy released from 
this chemical reaction is used to generate conforma- 
tional changes of the motor proteins that yield mechan- 
ical forces along cytoskelctal filaments. The dynamics 
of the resulting polymer network is controlled by active 
process on a range of time scales, including the poly- 
merization/deplymerization of the polar filaments, the 
force-generation form crosslinking motor proteins, and 
the load-dependent dynamics of these active crosslink- 
ers. 

Theoretical work has modeled the cytoskeleton via 
generic continuum hydrodynamics as an active liquid, 
where the effect of activity is incorporated via suitable 
modification of the hydrodynamic equations of equilib- 
rium liquid crystals [3H3 . The continuum theory has led 
to several predictions, including the onset of spontaneous 
deformation and flow in active films [7] , the formation 
of spiral and aster patterns reminiscent of those observed 
in in-vitro extracts of cytoskeletal filaments and motor 
proteins[3l |4j [8-13 , and activity-induced thinning and 
thickening in sheared active suspensions [Tlrtl7| . Vis- 
coelasticity has also been incorporated in the continuum 
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theory using the Maxwell model that modifies the re- 
sponse of the liquid by introducing a characteristic time 
scale controlling the crossover from fluid behavior at long 
times to elastic behavior at short times [5] . Given, how- 
ever, that the active liquid viscoelastic model cannot sup- 
port elastic stresses at long times, its direct relevance 
for the understanding of the crawling dynamics of the 
lamellipodium and of active contractions in living cells 
remains to be established. In addition, the active liq- 
uid model is inadequate to describe cross-linked contrac- 
tile systems, such as stress fibers (cross-linked bundles of 
actin filaments and myosin minifilaments that play a cru- 
cial role in controlling the ability of non-muscle animal 
cells to generate and resist forces) [TS] or muscle sar- 
comeres that often exhibit spontaneous oscillations |19j . 
Such oscillations require long-wavelength elastic restoring 
forces [213 [2TJ not accounted for in an active (even vis- 
coelastic) liquid. This suggests that the long-wavelength 
properties of stress fibers or sarcomers may be better de- 
scribed as those of an active elastic medium or active 
solid. Polarity is generally expected to also play an im- 
portant role in these systems indicating that a suitable 
continuum model maybe that of an active polar elastomer 
gel. 

Passive polymer gels are often classified on the basis of 
the nature of the crosslinking forces [22]. Chemical gels 
have strong cross-links bound by covalent bonds. These 
crosslinks have an essentially infinite lifetime on all exper- 
imentally relevant time scales and the gel behaves elasti- 
cally at long times, with a finite shear modulus. At short 
times, however, dissipation induced by internal frictional 
processes can result in "liquid-like" response, with the 
loss (viscous) component of the elastic moduli exceed- 
ing the storage (elastic) component. In physical gels, 
in contrast, the crosslinks are held together by weaker 
interactions (e.g., dipolar or ionic) and have finite life- 
times, ranging from minutes to a fraction of a second. 
This yields a broad spectrum of behavior, from strong 
physical gels, that are similar to chemical gels, to weak 
physical gels, with reversible links formed by temporary 
associations between chains. The latter are liquid at long 
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time and exhibit elasticity on short time scales. 

Similarly, active polymer gels also may or may not ex- 
hibit low frequency elasticity, depending on the nature 
of the crosslinkers. Cross-linked reconstituted actin net- 
works exhibit some of the properties of strong physical 
gels and display large active stiffening driven by molec- 
ular motors [23]. MacKintosh and Levine [MJ [35] and 
Liverpool et al. [26] showed that elastic networks with 
contractile forces induced by myosin II motors, described 
as static force dipoles, can account for both the large scale 
contractility and stiffening observed in experiments. In a 
recent paper Giinther and Kruse [20] also demonstrated 
that a continuum theory obtained by coarse graining a 
specific microscopic model of coupled sarcomeres does 
yield oscillatory states, as observed ubiquitously in these 
systems, provided the load-dependent on/off dynamics of 
motor proteins is included in the hydrodynamic model. 
Motor proteins are also directly involved in controlling 
mechanical oscillations and instabilities in cilia and flag- 
ella [2"TH2"5] and in the mitotic spindle during cell divi- 
sion [3U]. In all these cases the elastic nature of the net- 
work at low frequency is crucial to provide the restoring 
forces need to support oscillatory behavior, i.e., these sys- 
tems are best modeled as active solids, rather than active 
liquids. 

In this paper we formulate a generic continuum the- 
ory of isotropic cross-linked active gels that incorporates 
the on/off dynamics of crosslinking motor proteins. Fol- 
lowing MacKintosh and Levine (24] [25], we model the 
gel as a two-component system composed of an elastic 
network coupled frictionally to a permeating fluid. The 
details of the model are given in section II. The active 
forces arising from motor proteins are incorporated phe- 
nomenologically through an active contribution to the 
stress tensor of the elastic network and are controlled by 
the load-dependent on/off dynamics of the motors. In 
section III we examine the hydrodynamic modes of the 
active gel and show that, as stated in Ref. [5D], a large ac- 
tivity can change the sign of the effective compressional 
modulus, yielding a contractile instability. Spontaneous 
oscillations are obtained in the regime of weak activity 
where the compressional modulus is softened by bound 
motors, but remains positive. In section IV we consider 
the case of an overdamped gel relevant to muscle fibers 
and show that it can exhibit propagating waves and os- 
cillatory instabilities as parameters are varied. A phase 
diagram summarizing the behavior is given in Fig. 2. 
In Section V we describe the macroscopic homogeneous 
response of the active medium as probed in creep exper- 
iments and by macroscopic rheology measurements. The 
two-component gel model exhibits viscous response on 
short time scales and elastic response at long times [31] 
even in the absence of activity, when the time scale con- 
trolling the crossover between these two responses is set 
by the ratio of the viscosity and the compressional mod- 
ulus of the network. Activity renormalizes the time scale 
controlling this crossover. Finally, we conclude with a 
brief discussion. 



II. HYDRODYNAMICS OF ISOTROPIC 
ACTIVE GELS 

Hydrodynamics is a systematic method to study the 
behavior of extended systems on long times and length 
scales by focusing on the dynamics of conserved and 
broken symmetry fields. Here we use a phenomcno- 
logical symmetry-based approach to formulate a con- 
tinuum hydrodynamic description of a cross-linked gel 
(e.g., a network of actin filaments crosslinked by filamins 
or other "passive" linkers) under the influence of active 
forces exerted by clusters of crosslinking motor proteins 
(e.g., myosin II minifilaments) . We consider a three- 
dimensional isotropic polymer gel of mesh size £, vis- 
cously coupled to an incompressible permeating New- 
tonian fluid |31| . This two component model has been 
used previously to determine viscoelastic response of a 
filamentous isotropic network in solution [24, 25, 31, 32], 
and more recently to discuss mechanical response of a 
coupled network-solvent system when probed by an ac- 
tive agent [33]. At length scales larger than £ the de- 
formations of the polymer network can be described by 
isotropic elasticity in terms of a continuum displacement 
field, u(r, t) and an elastic free energy given by 

Fe = \j ( Xu *i + fyUijUij) , (1) 

with A and p, the usual bulk and shear Lame coefficients 
and Uij = ^(diUj +djUi) the strain tensor. The permeat- 
ing viscous fluid is characterized by a velocity field v(r, t) 
and the coupling between the network and the fluid is 
controlled by a friction per unit volume, T. The equa- 
tion of motion for the displacement field can be written 
as 

pa = -r(u - v) + V • a , (2) 

where p is the mass density of the network and <r is the 
stress tensor of the gel. The permeating fluid is described 
by the Navier-Stokes equation, 

p/v-r7V 2 v + VP = T(u- v) (3) 

where pf is the mass density of the fluid, rj the fluid 
shear viscosity, and P is the pressure. We have assumed 
a low Reynolds number regime for the fluid and omitted 
the convective term from the Navier-Stokes equation. It 
is also assumed that motor proteins do not exert any 
direct forces on the permeating fluid. As discussed else- 
where |31) , the friction T between the elastic network and 
the permeating fluid can be estimated by considering a 
polymer strand of length £ moving relative to the back- 
ground fluid at a velocity v. By equating the viscous 
force density ~ t/v/£ on the strand to the viscous fric- 
tion ~ Tv due to the permeating fluid, one obtains an 
estimate of the friction as T ~ v/S, 2 - The frictional drag 
per unit volume T is then determined by the force density 
required to drive a fluid of viscosity r\ through network 
pores of characteristic cross section £ 2 . 
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The stress tensor of the gel can be written as the sum 
of elastic, dissipative and active parts, 



a = <t + a + <j 



(4) 



The elastic contribution is given by af- — with 

<*y = ( A + y ) 5 « V • u + 2 M («« - |<5ij V ■ u) . (5) 

The dissipative component cr d is given by 

<rfj = Vb^ijV -u + 2t] s (iL i:j - -SijV ■ u^j , (6) 

where rjb and rj s are bulk and shear viscosities arising 
from internal friction in the gel. Changes in the density 
p of the network are slaved to changes in volume, thus 
5p = — poV-u, with /?o the mean mass density of the elas- 
tic network. In addition, we neglect here for simplicity 
energy fluctuations and assume that the fluid surround- 
ing the network serves as a heat bath and maintains the 
temperature constant. This approximation is not ade- 
quate to describe real muscle fibers that heat upon con- 
traction. 

The active contribution, <r a , to the stress tensor arises 
from the forces exerted by motor proteins bound to the 
filaments. We assume a total concentration c — C5 -\- c u 
of motor proteins in the gel, with ct, and c„ the concen- 
trations of bound and unbound motors, respectively. In 
an isotropic network the active contribution to the stress 
tensor can generically be written as [4], 



time T on that the protein spends attached to the filament, 
doing its working stroke, and the time r Q // that it spends 
detached from the filament, making its recovery stroke. 
Motor proteins are generally characterized by the value 
of the duty ratio, r = T on j(r on + r a ff). Myosins II with 
r ~ 0.05 [2j spend most of their time unbound, while 
two-headed kinesins have values of r close to unity and 
are classified as highly processive motors that remain at- 
tached to the filament for most of the duration of the cy- 
cle. The binding and unbinding rates are then estimated 
as k u ~ 1/t „ and kb ~ l/T tr. For individual myosins 
II, T on ~ 2ms and r ff ~ 40ms [2J, corresponding to 
k <C k u . During the working stroke and in the absence 
of external load, the protein moves along the filament at 
a speed v ~ Ap. The time r on ~ 1/vq ~ 1/A/z depends 
on motor activity, while r a ff is essentially independent 
of A/i. For myosins II this gives k u ~ Ap kb- 

Finally, we assume that the four-component active gel 
described by the set of coupled equations ([3]), ^ and 
Q is incompressible. This requires 



V • [(1 - (j> p ) v + (j> p u} = , 



(10) 



< = <% C(p> c b) &p 



(7) 



where <p p denotes the combined volume fraction of the 
polymer network with bound motors. We assume that 
the volume fraction of the network is very small, i.e. 
4> p « 1. In this case Eq. (10) reduces to the condition 
of incompressibility of the ambient fluid, V • v ~ 0. 

In the homogeneous steady state the network and fluid 
densities have constant values po and pf, respectively. 
The relative concentrations of bound and free motors are 
controlled by the binding/undinding rates and are given 
by 



where Ap is the change in chemical potential due to the 
hydrolysis of ATP and C(pi c b) is a scalar function with 
dimensions of number density describing the stress per 
unit change in chemical potential due to the action of 
active crosslinkers. 

To complete the hydrodynamic description we need 
equations describing the dynamics of bound and unbound 
motors. We assume unbound motors diffuse in the per- 
meating fluid, while bound motors are convected with the 
polymer network. Their dynamics is controlled by pre- 
scribed binding and unbinding rates, kt and k u according 
to first-order reaction kinetics. The resulting equations 
are 



d t c b + V • (c b u) 



-k u Cb + ki,c u 



(8) 



Cba ~ K + k Cm0 ' 



C U 



kh ~t" 



(lla) 

(lib) 



with c m o the total steady state concentration of motor 
proteins. In the following we are mainly interested in 
non-processive motors like myosins II that are mostly 
unbound on average, with Cbo << c„o- In this case we 
neglect the dynamics of free motors that essentially pro- 
vide a "motor reservoir" and assume that c u ~ c u q in 
([8]). In addition, we expand C(Pi c fc) to linear order in 
fluctuations of the network density and motor concen- 
tration from their equilibrium values, Sp = p — po and 
5c b =Cb- Q,o, as 



d t c n 



D\7 2 c v 



+ k u Cb kbC u 



(9) 



where D is the diffusion coefficient for free motors. The 
rates fc& and k u depend of course on the specific type of 
motor protein considered. Each motor protein undergoes 
a conformational transformation during a cycle fueled by 
a chemical reaction, generally the hydrolysis of ATP [5J . 
The total cycle duration is determined by the sum of the 



C(p,c fc ) = Co + Ci — 1- C2 — 

P0 CbO 



(12) 



The microscopic parameter Co is related to a stall force, 
but will not play a role in the following. The parameter 
Ci arises from spatial variations in the motor density. 
Both Ci and £2 are expected to be positive for contractile 
systems. 
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III. HYDRODYNAMIC MODES AND LINEAR 
STABILITY OF HOMOGENEOUS STATIONARY 
STATE 

In this section we consider the linear stability of the 
homogeneous stationary state, with u = v = 0, p = pa 
and Cf, = Cbo by examining the hydrodynamic modes of 
the incompressible gel. The fluid density pf is fixed due 
to the condition of incomprcssibility. Using (12 1 for the 



active parameter £, the linearized hydrodynamic equa- 
tions are given by 

Pqu - pV 2 u -(X + fi- Ci A/i)V(V • u) = T (v - u) 



+rj s V'u + rib 



V(V-u) + C 2 AmV0, (13a) 

PfV - r?V 2 v + VP = T (u — v) , (13b) 
+ V • ii= -k u <p , (13c) 

with <p = dcb/cbo and the condition V ■ v = 0. 

We now discuss the hydrodynamic modes of the three- 
component system described by Eqs. ( 13a|13c I obtained 
by neglecting fluctuations in free motors. We expand 
the fluctuations Sy a — (u, v, <fi) in Fourier components 
according to 



-«q-r 



Sy a (fi,t) 



(14) 



and look for solutions with time dependence of the form 
Sy a (q, t) ~ e~ tut 5y a (q). We also write u into its com- 
ponents transverse and longitudinal to q by letting u = 
qitL + ut, with q = q/q and q • ut = 0. Due to in- 
comprcssibility of the background fluid, v does not have 
any longitudinal component, and incompressibility allows 
us to eliminate the pressure P from (13b). In Fourier 
space, dropping for simplicity of notation the tilde on the 
Fourier components of the fluctuations, the equations for 
the longitudinal fluctuations are given by 



p UJ 



CiA/V - 

(— ilj + k v 



->(r + VlQ 2 )] u l = -iqC,2&.p,(t> 
(15a) 

p = ujqu L , (15b) 



where we have defined the longitudinal modulus of the 
gel as B — A + 2p and a longitudinal viscosity of the 
network as tjl = Vb + (4/3)ry s . The longitudinal part of 
the displacement couples to motor density, but not to 
the velocity of the permeating fluid in the incompressible 
limit considered here. Fluctuations in the longitudinal 
displacement are slaved to fluctuations in the network 
density, with dp — poiqu^. The longitudinal equations 
(15a I and (15b) can then also be rewritten as coupled 



equations for fluctuations in the network and bound mo- 
tor densities, 



[- PQ u 2 + (B- CiA/V 



(— ilu + k n 



:(r + VL q 2 )} S -P=( 2 Ap q 2 t., 
Po 

(16a) 

i = —lui— , (16b) 
Po 



Finally, the equations for the transverse components are 
given by 

[-pou 2 - ioj{T + Tjsq 2 ) + pq 2 ] u T = Tv , (17) 



(—Kdpf + r + rjq 2 ^) v = — iuiTut . 



(18) 



and are decoupled from the equations for the longitudinal 
modes. We therefore proceed to analyze the two groups 
separately. 



A. Longitudinal Modes 

In the incompressible limit considered here, the only 
role of the permeating fluid is to provide the frictional 
damping T. The longitudinal deformations of the poly- 
mer network do, however, couple to fluctuations in the 
bound motor density. It is instructive to first review the 
behavi or of a passive gel, as obtained by letting Ap = 
in Eq.( 15a). 



1. Passive gel. 

In the absence of motor proteins, longitudinal fluctu- 
ations in an incompressible gel are controlled by a single 
equation, given by 



{- 



■p u 



B q 2 - iu [r + W] } U L = 



(19) 



We stress that this equation also describes the behavior 
of fluctuations sin the network density, as 5p = iqp v,L. 
The hydrodynamic modes are the roots of the quadratic 
polynomial in curly brackets in Eq. (|19[) and are given by 



2po 



2\2 



±PoBq 2 



(20) 

The behavior is con trolled by the interplay of two length 
scales, £d = yJrjijT, the length scale over which intrinsic 
viscous dissipation within the network is comparable to 
dissipation due to friction with the permeating fluid, and 
£r = 2^/pqB/T controlling the ratio of elastic restoring 
forces in the network to viscous drag from the permeat- 
ing fluid. The length scale It has been introduced before 
by Mackintosh and Levine [24l [25]. At small wavevec- 
tor (q <C ^p -1 ) the dispersion relations are always imag- 
inary, corresponding to relaxational or diffusive modes, 
and take the form 



"r 






— h 

Po 


(- 

\Po 





B 2 
- l Y q 



0(q 4 



+ 0(q 4 ), (21) 
(22) 



where the superscript is used to denote the passive 
gel limit. The mode ui L r is non- hydrodynamic and de- 
scribes the relative motion of the polymer network and 
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the permeating fluid. The mode w° describes the dif- 
fusive relaxation of network density fluctuations. In the 
two-fluid incompressible gel model considered here there 
are no propagating longitudinal sound waves |3f) and the 
network density Sp relaxes diffusively, while the solvent 
density pj remains fixed. The limit q <C ^p 1 holds if 
< £d- On the other hand, when £p > the modes are 



relaxational as given in Eqs. (21 ) for q <C t T , but there 
is an intermediate regime of 1^ < q < £7* where the gel 
can support propagating sound-like density waves. Prop- 
agating density waves exist if the argument of the square 
root on the right hand side of Eq. (20) is positive, i.e., 
ApoBq 2 > [r + rji^q 2 ] 2 . It is convenient to scale lengths 
by £d, with q = q^d- The condition for the existence of 
propagating waves can then be written as 

(l + eq 2 ) 2 



B* > 



4q 2 



(23) 



where B* = and e = til/t). The propagating waves 

are controlled by the interplay of inertia and elasticity 
and decay on time scales of the order of the relaxation 
time Tr = po/T, which is set by the frictional damping 
from the solvent. The equality sign in Eq. (23) defines 
the critical line shown in Fig. [I] separating the region 
of diffusive density relaxation from the region where the 
system supports propagating sound-like waves. No prop- 




Diffusive Modes 
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FIG. 1: The critical line B*{q) given in Eq. (23) for e = 1 sep- 
arating the region of parameters where density fluctuations in 
a passive incompressible gel relax diffusively, from the region 
where the system supports propagating density waves. In the 
chosen dimensionelss units, the same line also describes the 
boundary £>* (q) obtained for the case k u — > oo and boundary 
B'a(q) obtained for the case k u = 0. 

agating waves exist for B < T 2 ^ 2 /po, corresponding to 
the minimum of the curve in Fig. [T] We stress that the 
modes are always diffusive at the longest wavelengths, 



when q — > 0. These finite wavevector sound-like waves 
persist down to very small wavevector in the limit of van- 
ishing friction T with the surrounding fluid. This is seen 
by setting T = first, followed by the small wavevector 
approximation. The dispersion relations then take the 
form 



±Q\ n Ty- 
Po 2p 



(24) 



These are indeed sound waves propagating at the longi- 
tudinal sound speed ~ y^B/po. 



2. Neglecting bound motor fluctuations (k u — > oo). 

We now proceed to incorporate the effect of motor 
proteins. We first consider the case of stationary bound 
motors. This can be obtained in two ways, either by 
letting k u — > co, which corresponds to neglecting bound 
motor fluctuations, or by letting k u — 0, which corre- 
sponds to neglecting the motor on/off dynamics. In both 
limits motor activity can yield a contractile instability of 
the system, but no spontaneous oscillations, as pointed 
out in Ref. |2"TJ] . 

When k u — > oo, then = and the concentration of 
bound motors is constant, cj = q,o- We then obtain a 
single decoupled equation for fluctuations in the longitu- 
dinal displacement (or equivalent, in the network density 
8p) of the form 



{- 



'PqU 



(B-^I\p)q 2 ~itu[T + r, L q 2 }}u L . (25) 



In this limit the only effect of motor activity is a con- 
tractile reduction of the compressional modulus, which 
is given by 



B a = B- CiAp 



(26) 



The hydrodynamic modes are identical to those de- 
scribed in the previous subsection, with the replacement 
B —> B a . If B a < the imaginary part of the mode 
luI changes sign, signaling a contractile instability of the 
system driven by motor activity. When B a > 0, the 
modes can be real at finite wavevector, corresponding to 
propagating waves. The condition for the existence of 
propagating waves is precisely as given in Eq. ( 23 ) for 
the passive gel, with the replacement B — > B a . A plot of 
B* = B a p /(T 2 t; 2 ) as a function of q is that identical to 
that shown in Fig. [T] for the passive case. We stress that 
the existence of these propagating density waves is not 
a consequence of activity. There is in fact a maximum 
value of activity, given by £iAp c = B — T£ 2 /ty, and cor- 
responding to the minimum of the curve plotted in Fig. [T] 
above which there are no propagating modes. In addi- 
tion, since £r ~ \J B a decreases with increasing activity 
Ap, the range of wavevectors where propagating waves 
exist for a fixed B* decreases with increasing activity and 
is given by Aq = 2{B* a (B* a - I)} 1 / 4 . 
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3. Neglecting bound motor dynamics (k u — 0). 

In this case bound motors remain bound at all times 
and bound motor fluctuations are slaved to network den- 
sity fluctuations, with 4> — Sp/ pq. The relaxation of lon- 
gitudinal fluctuations is described by 

{-Pocj 2 + [B - (Ci + C 2 )A/i] q 2 - ilo [r + n L q 2 ] } u L = . 

(27) 

and the only effect of static bound motors is a further 
downward renormalization of the elastic modulus, which 
is now given by 



B' a = B - (Ci + C 2 )A M 



(28) 



The modes are again formally identical to those obtained 
for the passive gel, but with B — > B' a . The gel exhibits 
a contractile instability for B' a < and finite wavevector 
progating density waves for B' a > 0. Note that the limit 
where all motors are bound can be obtained for instance 
after full hydrolysis of ATP to ADP. Myosin has a high 
affinity to actin, hence in a pure ADP environment it will 
act as a "permanent" bound crosslinker [33]. In this case, 
however, there will also be no reduction of the elastic 
modulus due to activity, hence no contractile instability. 
In fact muscles become rigid as ATP runs out, which is 
one of the causes of rigor mortis. 



4- Including bound motors dynamics (finite k u ). 

We now incorporate the dynamics of the bound mo- 
tors and consider the hydrodynamic modes of the the two 



coupled equations (I6a) and (16b). These yield a cubic 



eigenvalue equation, given by 

ip uj 3 - lu 2 (Y + k u p + riLq 2 ) - iw|fc u (r + r) L q 2 ) 
+ [B- (Ci + Ca) q 2 } + k u B a q 2 = . (29) 

The behavior is now controlled by the competition of 
two time scale, the network relaxation time rp = Po/r 
and the time scale r on = fc,^ 1 characterizing the motors 
on/off dynamics. Solving perturbatively for small wave 
numbers q, the three modes are given by 



= -i^q + 0{q ) , 



iq 



C2A/X 
1 - Kpo 



+ 0(q i ), 



(30) 



(31) 



wavevector. Again, it changes sign when Ap > -B/Ci, 
corresponding to a contractile instability of the gel that 
occurs when the active stresses exceed the elastic restor- 
ing forces from the passive elements of the polymer net- 
work. The other two modes are non-hydrodynamic and 
always stable at long wavelengths. The mode with re- 
laxation rate uj b describes the decay of fluctuations in 
the density of bound motors. The mode with relaxation 
rate ujl,v describes the damping of the network due to 
its motion with respect to the permeating fluid. Even 
when the on/off dynamics of the bound motors is taken 
into account, no spontaneous oscillations are generated 
by motor activity in the long wavelength limit. Oscilla- 
tory solutions do, however, occur at finite wavevector, as 
described below. We note that, although the modes al- 
ways remain stable, the coupling to motor activity can 
yield a change in sign of the 6(jq 2 ) damping in LOb and 
(jJl,t- This effective "negative viscosity" due to motors 
occurs when the time scale of the motor on/off dynam- 
ics is fast compared to the frictional relaxation of the 
network, i.e., for Tr > r on . This "negative friction" ef- 
fect of motors will become important below and was also 
discussed in Prost et al [55] . 

As in the passive case, the dispersion relations of the 
hydrodynamic modes of our model viscoelastic gel de- 
pend on the order in which the limits Y — > and q — > 
are taken. Above we considered the small q limit for fixed 
r. If in contrast we take r — > first, followed by q — s- 
we obtain propagating modes (for B a > 0). The mode 
LUb describing relaxation of bound motor fluctuations is 
qualitatively unchanged and takes the form 



Lu b = -ik u - n 



: C 2 Ap 
k u po 



0(q 4 ) 



(33) 



The two modes lot, and WL,r describing the dynamics of 
network density fluctuations are replaced by two propa- 
gating modes (for B a > 0), with dispersion relation 



WL,± 




C 2 A M 
ku 



(34) 



In contrast to the case of a passive gel or a gel with static 
bound motors, these oscillatory density waves can now 
become unstable when the (negative) viscosity induced 
by the motors overcome the internal viscous dissipation 
of the network, i.e., for (2A/it „ > r/L- Above the critical 
value of activity defined by the vanishing of the damping 
in Eq. ( 34 1 , the propagating waves become unstable and 



the uniform state is presumably replaced by a state that 
supports spontaneous oscillations. 



r 2 f B a C2A/X T)l\ ^ f 3N B - Transverse Modes 



Po V r r - k u p p ) 

(32) The transverse equations ( [l7| and ( 18 ) do not couple 

The mode lot, describing the relative mass diffusion of to motor dynamics. They yield a cubic eigenvalue equa- 

network and solvent in the gel is unchanged at small tion. There are therefore three transverse modes in the 
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system. Of these two are propagating shear waves, with 
dispersion relation for small q given by 



w(q) = ±q 




r/ + Vs + 




0(q 3 ) , (35) 



with p g = po + pf the mass density of the gel. The 
third transverse mode is a non-hydrodynamic mode with 
a finite decay rate at q = 0. It describes the relative 
motion of the polymer network and the permeating fluid. 
The dispersion relation is given by 



w(q) 



iVpg_ 
PoPf 




Transverse fluctuations always decay and to linear or- 
der do not destabilize the stationary homogeneous state. 
Finally, if B and p are comparable, the speed of propaga- 
tion of the transverse waves given in Eq. ( 35 ) is generally 



much smaller than that of the longitudinal waves given 
in Eq. (24), since pf » po- 



IV. OVERDAMPED DYNAMICS OF THE 
POLYMER NETWORK : CONNECTION TO 
MUSCLE SARCOMERES 



In the overdamped limit of large friction T 
tial term in Eq 



(13a 



the iner- 

is negligible and the relaxational 
dynamics of the fiber density is controlled by the viscous 
coupling to the permeating fluid. This is the limit that is 
relevant to most biological systems, such as muscle sar- 
comeres. We show here that in this limit the on/off dy- 
namics of bound motor yields an effective inertia that re- 
sults in spontaneous oscillations even in this overdamped 
limit. 

The approximation of neglecting the inertial terms can 
be quantified as follows. The inertial term in Eq. ( 13a 



can be neglected relative to the frictional damping from 
the fluid provided pquj 2 << Tu> or u> << T/po ~ 
rj/(f po)i which is simply the condition of low Reynolds 
number for an object of typical size £ moving in a medium 
of kinematic viscosity ry/po at a typical speed ~ 
A sarcomere of typical rest length £ ~ 2.5 pm [I], 
moves in an ambient viscous medium of viscosity rj ~ 
10 pN s pm^ 1 . The mass density po of a sarcomere is 
approximately 10 3 kg m~ 3 f35]. Inertial effects can be 
neglected if the velocity of a sarcomere unit, typically of 
order 10 pm s" 1 , is small compared to ??/£po- From the 
known values of sarcomere parameters, as quoted above, 
v/^Po ~ 10 ms , which is three orders of magnitude 
higher than the typical velocity of a sarcomere. Hence 
the ignoring of the inertial forces is justified. 

A sarcomere chain can be described as a one dimen- 
sional elastic system in terms of a displacement field 
u(z,t), with z the coordinate along the sarcomere's 
length. In the overdamped limit the equation for the 



displacement field and the deviation of the fraction of 
bound motor from the steady state value are given by 



(r 

d t <t> 



VL d 2 ) d t u = B a d 2 z u + ( 2 Apd z 

= -0,[(i + 0)ft«]-M- 



(37) 
(38) 



We note that in the overdamped limit discussed in this 
section our model is formally similar to the model intro- 
duced by Murray and Oster [31| to describe the role of 
the mechanochemistry of the cytogel in epithelium move- 
ments (albeit with calcium dynamics taking the place 
of motor dynamics), but with one important difference: 
here we consider a gel frictionally coupled to a perme- 
ating fluid, while Refs. [23 [35] consider a gel elastically 
coupled to a substrate. As shown below, both models 
yield oscillations and traveling waves. 

When linearized by approximating the convective term 
on the right hand side of Eq. (38) as ~ —d z dtu, these 



equations arc identical to those derived by Gunthcr- 
Kruse [20] from a microscopic model of muscle sarcom- 
eres. Here we show that the same equations can be ob- 
tained by a purely phenomenological approach that in- 
cludes both the dissipation due to the coupling to the 
permeating fluid and the on/off motor dynamics. We 
also note that the bound motor fraction can be eliminated 
from the linearized equations by transforming them into a 
single differential equation for the displacement. Solving 
the linearized form of Eq. ( 38 ) for <f> with <p(z, t = 0) = 0, 



substituting in Eq (|37j) and differentiating with respect 
to time, we obtain a single differential equation for the 
displacement u(z,t), albeit second order in time, given 
by 



(r - rj L d 2 ) d 2 u + [r - i lL dl - Va d 2 z ] d t U 



where 



Va=T on [B- (Cl+&)A/i] 



B a d 2 z u 
(39) 

(40) 



It is clear from Eq. ( 39 ) that the effect of motor on/off dy- 
namics is to provide an "inertial" contribution to the dy- 
namics of the network. On length scales large compared 
to £d we can neglect the internal dissipation intrinsic to 
the network proportional to the viscosity rji, compared to 
the friction T with the permeating fluid. Eq. (39) then 
simplifies to 



T on Td 2 u +[T- Vad 2 z ] d t u = B a d 2 z u 



(41) 



In this limit Eq. (41) describing deformations of the ac- 
tive network is formally identically to Eq. ( 19 ) for the 
passive gel, with t ov T playing the role of a mass den- 
sity, and a viscosity rj a and an elastic modulus B a , both 
renormalized by activity. The effective viscosity and the 
elastic modulus can change sign at high activities, yield- 
ing instabilities. 

First we consider the hydrodynamic modes of the sys- 
tems described by the linearized form of Eqs. (37) and 



(38) or by Eq. (39). These are given by the solutions of 
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the eigenvalue equation, given by 

w 2 (r + r) L q 2 ) + iojk u [r + (r) L + Tja) q 2 ] - k u B a q 2 = 
The general solutions of the eigenvalue equation are 



(42) 



' 2(T + VL q 2 



^- [T + ( VL + Va ) q 2 ] 2 + ^^(r + VL q 2 ))m 



For small wavevector {q — ¥ 0) we obtain two modes, 



w 6 = -lk u + i - q 

B a 2 



(44) 
(45) 



describing motor and network density relaxation, respec- 
tively. Again, the system exhibit a contractile instability 
when B a < 0, but there are no oscillatory waves in the 
long wavelength limit. 

Propagating wave solutions exist if the argument of 
the square root on the right hand side of Eq. ( 43 1 is pos- 
itive. The active viscosity can be written as rj a = (B a — 
C2A/i)/fc„, hence it depends on the renormalized elastic 
modulus B a . If we choose to treat B a — B a / (T^k u ) and 
C2 = C2 A/^r^jfcn) as independent parameters the con- 
dition for existence of propagating waves can be written 
as B-(q)<B a <B+(q), with 



1 



i + <f + C 2 g 2 ±2V? 2 (i 



q 2 + (2q 2 ) 



(46) 

where we assumed t\l ~ T£ 2 . Propagating waves then ex- 
ist in a band in the {B ai q) plane, as shown in Fig. [2] The 

width of the band is AB a = A\J 1 + C2 + l/<7 2 - It van- 
ishes at small wavevectors and goes to the constant value 
4(1 + C2) 1 ^ 2 at large wavevectors. In contrast to the prop- 
agating density waves obtained in a damped passive gel, 
the oscillatory behavior results here from motor activity 
and the range of parameter where it exists grows with 
the time r on that characterizes motor dynamics. Since 
Ton ~ 1/A/x to leading order £2 is independent of activ- 
ity for small activity. In addition, the propagating waves 
are unstable when the imaginary part of the eigenvalues 
given by Eq. (43) is positive. This corresponds to 



B a < C2 - 



1 



(47) 



and defines a region where the overdamped active gel 
exhibits an oscillatory instability. We expect that when 
nonlinear terms are included in the equations, the gel will 
exhibit spontaneous oscillations in this region of parame- 
ters. The transition from diffusive to oscillatory behavior 
is controlled by the interplay between r OTl and the char- 
acteristic time Td ~ ^T/B for the diffusive relaxation of 



B,r 5 




Diffusive 



Propagating Waves 



Diffusive^ 



Oscillatory Instability 



Contractile Instability 



FIG. 2: A phase diagram for the overdamped active gel. The 
vertical axis is B a = B a /(r^k u ) and the horizontal axis is 
q£,d- The boundaries separating the regions of diffusive relax- 
ation of network density fluctuations from the region where 
traveling waves exist are given by Eq. ( 46 1 , plotted here for 



£2 = 2. Below the horizontal line B a = 0, the system exhibits 
a contractile instability. 



a network fluctuation of size If r OJl 3> Td the on/off 
motor dynamics provides an "inertial drag" to the net- 
work that opposes the elastic restoring forces, yielding 
propagating waves. Alternatively, the result can be un- 
derstood in terms of two length scales in the problem, ^ 
and lb ~ yf B / (k u T). If ^ > li then density relaxation 
is always diffusive in the range of wavevectors {q^d <C 1) 
described by the present theory. If in contrast lb > £<j 
the network supports propagating density waves in the 
wavevector range 17 < q < £7 . 



V. LINEAR RESPONSE 

A. Dynamic Compressional Moduli 

In this section we characterize the macroscopic ho- 
mogeneous viscoelastic response of the active gel in fre- 
quency space in terms of the dynamical compressional 
modulus. To describe a traditional compressional exper- 
iment, we consider a slab the three-fluid active gel model 
with only longitudinal degrees of freedom, held between 
two plates at z — and z — L and unbounded in the 
other two directions. We imagine applying a harmonic 
compressive strain at one end, where u{z = L) = u^e^ luJt , 
while holding the other end fixed, i.e., u(z = 0) = 0. In 
general, both the cases of an oscillating boundary that is 
permeable or impermeable to the permeating fluid are ex- 
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perimentally relevant. To implement a calculation that 
allow to treat both cases one needs to include a finite 
compressibility so that the longitudinal elasticity equa- 
tions couple to the fluid velocity v. Here we limit our- 
selves to a permeable boundary and impose no bound- 
ary conditions on v. With these boundary conditions we 
calculate the stress a(z = L) required at the oscillating 
boundary and define the complex compressional modu- 
lus B expt (uj) measured in experiments as the ratio of the 
stress to the applied compressional strain, Uq/L. We will 
see below that at low frequency we recover the complex 
bulk compressional modulus, B(uj) = B — ilotil, obtained 
assuming an affine compression over the entire sample. 

First we analyze for comparison the case of the passive 
gel with inertia and damping. The elastic response is 
governed by the equation 



p dfu + Td t u = Bd 2 z u + r^ L d t d 2 z u 



(48) 



We assume a solution of the form u{z,t) = f(z)e~' luJt , 
where f(z) = fie XiZ , yielding a characteristic equation 
for the eigenvalues A, 



A 



ui 2 po + uoY 

W) 



(49) 



Boundary conditions, /(0) = and f(L) = uq lead to 
the solution, 



. sinh(Az) 
f(z) = u — 



(50) 



sinh (XL) 

The complex dynamic compressional modulus is then 
given by B exp (co) = j^B(u) (|Q _ which gives 

B expt {uj) = B(uj)XL coth (XL) (51) 
The eigenvalue can be written as 



X 2 L 2 



ujL 


2 

+ i 


L 


_v s (uj)_ 




[6(u>)\ 



(52) 



where we have defined the frequency dependent sound 
speed, v s (uj) = yj 'B(ui) / po, and the penetration depth 
S(uj) = *J B(uj) /ojT which controls the penetration of rar- 
efaction/compression waves of frequency w [39] . At low 
frequency, where \XL\ <C 1, we recover B exp t(oj) — } B(uj) : 
provided ojL/v s (uj) <C 1 and L <C S(uj). The first condi- 
tion means that the frequency of applied oscillations is 
small compared to the frequency of sound wave propaga- 
tion across the entire sample. When this is not satisfied 
there is an appreciable time lag between the imposed 
deformation at one end of the sample and the deforma- 
tions realized at other material points across the sam- 
ple, resulting in nonuniform strain and preventing the 
experimental determination of a macroscopic compres- 
sional modulus. The second condition demands that the 
boundary compressional waves fully penetrated the sam- 
ple, which is again necessary to achieve a uniform com- 
pressional strain. For a similar discussion of shear rheo- 
logical experiments see Appendix C of Ref. [3T] . Finally, 



Z& pt M/rf/fc„ 

1000 
10 
0.1 
0.001 



ft "(a>) 



0.001 0.1 



10 1000 



FIG. 3: Storage (B' expt (uj)) and loss (B" xpt (u))) moduli, for 
B a = 1.15, C 2 = 1.1 and £ d /L = 0.5. 



the compressional modulus to second order in frequency 
as measured in a macroscopic experiment is given by 

B expt (uj) = B-u J 2 p L 2 /3-iu J (r lL +rL 2 /3)+O(u J 3 ) (53) 

We now turn to the compressional response of an active 
gel. In this case we ignore the inertial contributions rel- 
ative to the damping from the permeating fluid and look 
for solutions of the linearized version of Eqs. (|37| and ( 38 ) 
of the form u(z = L) = u e~ 4W * and 4>(z,t) 
with f(z) = fie x ' z and g(z) = 
are given by 



g(z)e-^, 
The eigenvalues 



X 2 L 2 



L 


2 r 


_Sa(u)_ 





1 



IUJ( 2 AP/Bg(uj) 

— IL0 + ki, 



(54) 



where S a (u) = yjB a (u)/wT and B a (uj) = B a - iur\ L - 
Using, —iLuf'(z) = (—VjJ + k u )g(z) and the boundary 
conditions on f(z) and proceeding as in the passive case, 
we obtain 



Bexpti^) - 



B a (uj) + C2A/1: 



lUJT on 



lU!T n 



XL coth (XL) (55) 



The real and imaginary parts of B^ xpt (uS) = B' expt (ui) — 
iB" xpt (w) representing the storage and loss moduli, re- 
spectively, are shown in Fig.(|3| for generic values of pa- 
rameters. The storage or elastic modulus has a frequency 
independent plateau at frequencies lower than the mo- 
tor's unbinding rate, indicating that the system behaves 
like an elastic gel in this region. The linear frequency de- 
pendence of the loss modulus is the hallmark of a dissipa- 
tive gel. At low frequency B' expt > B" xpt and the system 
behaves elastically, while at high frequency B" xpt > B' expt 
and the response is dominated by viscous losses. This 
response is reminiscent of the Kelvin- Voigt model of vis- 
coelasticity. Finally, at low frequency the compressional 
modulus is given by 

/ YL 2 

Bexpti^) = B a - UJ 2 T 2 n ( 2 A[l - iLJ ll] L + — T on ( 2 Ap 



-O(to 3 ) 



(56) 
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we obtain 

B a expt {u) ~ (B a - C 2 A^ - iur, L ) (L/U) coth (L/&) 



(57) 



whereas, at high frequencies since A ~ y/T/r)L — viscoelastic equation [43] . In this case the creep following 

application of a sudden load at t = 0, a{t) = <To0(t) has 
the familiar form 

«(«)= 8 :°!i n (i-e-^ L2 ^), 

"f]L^ +11/ 

where = t]l/ B a is the Kelvin- Voigt relaxation time. 

For finite values of r on , the creep response is controlled 
by the interplay of the two times scales tb and r on . We 
assume B a > 0, corresponding to weak activity. When 
B a < the system exhibits a contractile instability and 
the strain becomes arbitrarily large at long times for any 
applied cr{t). The evolution of the strain in response to 
an applied stress is then controlled by the two eigenvalues 
of Eq. ([59| for m = 1, given by 



B. Creep 

Here we study the macroscopic behavior of our ac- 
tive elastic medium by considering the creep response, 
i.e., the time evolution of the average strain e{t) = 
1/L Jq dz d z u in response to a homogeneous external 
stress, cr{t). In particular we are interested in character- 
izing the load and recovery creep of the material follow- 
ing the sudden application and removal, respectively, of 
a constant stress. Both responses are measured experi- 
mentally in cells [40 -42 . 

Consider a muscle fiber of length L with free bound- 
ary conditions at the ends z = and z = L, i.e. 
d z u{z = 0, L) — 0, and no fluctuation in motor den- 
sities being imposed at the ends. Hence one assumes 
normal mode expansions for u and (f> to be of the 
form u(z,t) = J2m=o u m{t) cos {rhz), and (j>(z,t) = 
Z)m=i ^rn{t) sin {rhz), where rh = mir/L. 

Neglecting nonlinearities, the evolution of the normal 
modes u m {t) and 4> m {t) in the material in response to a 
small external stress a(t) is governed by the equations 



A± = -7± \ h 



TB 



J{L)r on 

where time is measured in units of tb, and 

C 2 A/j | t b | TL 2 
B n 



(61) 



2 7 = 1 - 



/J(L), 



J(L) 



L 



(58a) 
(58b) 



The linear creep response of the active gel can then be 
classified as follows: 



I. 7 > , 7 > Jjj^r 
, 2 



stable monotonic behavior 



(r + rj L m 2 )u m (t) + B a m 2 u m (t) - C, 2 ^mcj) m (t) = f m (t) 
<j) m (t) = mu m {t) - k u <p m (t) . 

With, f m {t) = ^jt^ /(f dz sin {rhz) . 

Eliminating the fluctuations <t> m {t) in the density of H. 7 > , 7 < 7(7^ : stable oscillatory behavior 
bound motor, we obtain an effective equation for u m {t), 

given by m~ 7 — : sustained oscillations 

r on (r + VL rh 2 )u m {t) + [rh 2 {r on {B a - C 2 A^) + Vl) + T] u m {tf V - 7 < , 7 2 > 7(fc : unstable oscillatory growth 

+B a m u m {t) = Tonfm + fm ,y. 7 < , 7 2 < jrzyp — = unstable monotonic growth 

where r on = k~ l . The decay rates of the individual 
modes are 



2-y{m) 



B a - C 2 A M 



riL 1 



L 2 



(60) 



Eq. (60) shows that 7(771) is an increasing function of 
m, hence the higher modes decay at a faster rate. For 
simplicity we then consider only the first mode, m = 1. 
Thus we approximate the averaged strain developed in 
the material as e{t) ~ —2u\-k/L. Also note that neglect- 
ing viscous coupling to the fluid V amounts to considering 
the limit of the fastest mode m — > 00. 

In the limit r „ — > 0, when motors are unbound at 



all times, Eq. (59) reduces to the familiar Kelvin- Voigt 



The behavior is summarized in the phase diagram of 
Fig. [^displaying the various regions in the (r on , C2) plane 
for fixed J{L). We note that when t oji /t b is increased 
for fixed C2A/J/ B a > 1, the material eventually becomes 
unstable to stretching. Figures 5(a) and 5(b) show the 
time evolution of the strain in response to a step stress 
of height <T and duration T, a{t) = a [Q{t) - 9(t - T)], 
with initial condition e(0) = 0. The response in region 
I of stable monotonic decay is similar to conventional 
Kelvin- Voigt response. In region II of stable oscillatory 
decay the interplay of the two time scales T on and t b 
yields the possibility of a strain overshoot. For finite r on 
we also need to specify an additional initial condition 
determined by the initial distribution of bound motors, 
since e(0) = 0(0) - <t>{0)/r on . Fig. 
response in region III of sustained oscillations. 



5(b) displays the 
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(b) 



VI. DISCUSSION AND CONCLUSIONS 

We have presented a generic continuum theory of ac- 
tive gels, modeled as a viscoelastic solid with bound mo- 
tor proteins that induce active stresses in the medium. 
In the limit where the inertia of the network is neglected 
and the equations are specialized to one dimension, the 
model is equivalent to that proposed by Giinthcr and 
Kruse 20 by coarse-graining of a specific mechanical 
model of coupled muscle sarcomeres. For large values 
of the motor activity as measured by the rate of ATP 
consumption, A/x, the contractile action of bound mo- 
tors yields a diffusive (contractile) instability of the gel. 
This result has been obtained earlier in models of mus- 
cle sarcomeres [20] and actin bundles [44]. Here we show 
that it is a generic property of active elastic media. For 
smaller values of motor activity the interplay of solid elas- 
ticity and the binding/unbinding dynamics of the motor 
proteins yields propagating waves and eventually oscil- 
latory instabilities in the linear theory. Both stable and 
unstable oscillatory modes are obtained even in the case 
of an overdamped gel, as relevant to muscle fibers. We 
show that the finite time scale of motor on/off dynamics 
yields an effective "inertial" contribution to the dynam- 
ics of the elastic medium controlled by the time T on that 



motors spend bound to filaments (see Eq. (39)). One of 



the new results of the paper is the phase diagram dis- 
played in Fig. [4] for the macroscopic response of the sys- 
tem to external stresses. In the linear model sustained 
oscillations are only obtained for special parameter val- 
ues corresponding to a line in the (t „,^2) phase dia- 
gram. It is expected that nonlinearities neglected in the 
present work will have a stabilizing effect and replace the 



FIG. 5: (color online) Strain e(t) in response to a step-stress 
a(t) (dashed line, red online) with T = 10. Strain and stress 
are measured in units of ao and ao/B a , respectively and time 
is in unites of tb- The various curves correspond to differ- 
ent values of r on , spanning the regimes described above. Top 
frame: r on — (dotted line, green online), corresponding to 
a passive gel with Kelvin- Voigt response; t oii /tb = 0.2 (solid 
line, black online), corresponding to region I of monotonic sta- 
ble response; t „/tb = 0.5 (dashed line, purple online), and 
Tou/tb = 1 (dashed-dotted line, blue online), corresponding 
to region II of oscillatory stable response. All curves are for 
(,2/B a — 1. Bottom frame: r on /rs = 0.5 and C,2/B a — 3, 
corresponding to region III of sustained oscillations. 



unstable oscillatory response with stable self-sustained 
oscillations. The model considered is relevant for the 
description of motor-induced spontaneous oscillations in 
muscle sarcomeres and other active elastic media, and 
may provide a useful framework for the understanding of 
lamcllipodium crawling. 

We plan to extend this work in various ways. First, 
an analysis of the effect on nonlinearities is needed. Two 
classes of nonlinear terms are important in our model of 
an active gel. The first is provided by nonlinear con- 
vective terms in the equation describing the dynamics of 
bound motors, as shown in Eq. (38 1, and also includ- 
ing dependence of the unbinding rate k u on the elastic 



strain d z u developed in the gel. These are the simplest 
continuum manifestation of the highly nonlinear load de- 
pendence of the microscopic motor unbinding rate, which 
in turn plays an important role in controlling the motor- 
induced negative friction induced by the cooperative ac- 
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tion of motor proteins on biological systems elastically 
coupled to their environment [301 E3 ES] • A second class 
of nonlinearities arise from higher order terms in the ex- 
pansion of the active parameter £ given in Eq. 0. A pre- 
liminary estimate of the effect of these terms suggest that 
they stabilize the oscillatory growing modes and yield 
stable sustained oscillations. 

In the liquid state of an active system the polarity of 
actin filaments plays an important role. The coupling of 
polarity and flow has been shown to yield spontaneous 
flow [BJ, banded states of inhomogeneous concentration, 
and oscillatory states [7] . It is similarly expected that the 
coupling of polarity and elasticity will yield new phenom- 
ena in active solids, including spontaneous deformations 
and oscillations. To incorporate the effect of polarity we 
have begun to consider the properties of an active polar 



elastomer, where the orientational order can be induced 
either by elongated passive crosslinkers [46 or by the 
myosin minifilaments themselves. In addition, the lat- 
ter exert active force dipoles on the medium that induce 
active stresses coupled to the orientational order. A de- 
tailed discussion and analysis of such active elastomers is 
left for a future publication. 
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